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approximation to the solution of variable-coefficient 
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forces. When compared to the form of the exact solutions, 
this approximation is in terms of "elementary" functions. 
Two specific forms of attrition-rate coefficients are 
considered, allowing for different maximum effective ranges 
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be used to model a constant-speed attack against a static 
defensive position. It is shown that for these 
attrition-rate coefficients, the Liouville-Green 
approximation is not consistently reliable for predicting 
force levels, and yields exact results only under certain 
restrictive conditions. ' Furthermore it was found that 
methodology is not presently available to accurately predict 
from Liouville's normal form the error which will be 
incurred by invoking the approximation in a specific 
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I. 



INTRODUCTION 



although it is well recognized that combat is a complex 
random process, a model that provides "closed-form" 
solutions can often be an aid to the military analyst by 
explictly portraying interrelationships between combat 
parameters. Deterministic differential equations are 
frequently used to model the combat process, and vary 
greatly in their detail and level of complexity. The 
military analyst must be able to recognize significant 
variables and their relationships. An idealized model is 
useful for giving a "first cut" portrayal of such 
relationships. In particular, one can often obtain an 
analytical solution to such models. An analytical solution 
that can be expressed in terms of "elementary" functions can 
often provide insights into the combat process and identify 
factors affecting battle outcome. 

The use of differential equations to model warfare was 
pioneered by F. 0. Lanchester in 1914 [Ref 8]. This original 
work was quite idealized and assumed that the fire 
effectivenesses were constant throughout the battle. Many 
extensions to Lanchester' s classic formulation have 
subsequently been developed, to include such factors as 
mobility of forces and fire effectiveness that varies during 
the battle. Taylor and Brown [Ref 13] have developed a 
mathematical theory for solving Lanchester-type equations 
with these temporal variations in fire effectiveness. The 
fire effectiveness of forces is represented by what are now 
called Lanchester attrition- rate coefficients. Taylor and 
Brown give results for two specific forms of the 
coefficients: (1) effectiveness of each side's fire 
proportional to a power of time, and (2) effectiveness of 
each side's fire linear with time but a nonconstant ratio of 
attrition- rate coefficients. 
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Taylor and Brown present accurate numerical solutions 
to equations with the above attrition-rate coefficients. 
These solutions are, in general, stated in terms of infinite 
series. The solution method, however, is rather complicated 
and involved the development of new mathematical functions 
alien to the military analyst. Exact results are certainly 
beneficial, but the analyst also needs a solution in terms 
of "familiar" functions to permit parametric analysis and 
simplify recognition of significant variables. To provide 
such a solution form, Taylor [Ref 14] has suggested the use 
of the so-called Liouville-Green approximation, which can be 
expressed in terms of relatively simple functions. For an 
approximation to be useful in predicting force levels, 
however, it must be reasonably "close" to the exact 
solution . 

This paper provides a numerical evaluation of the 
Liouville-green approximation for the two forms of 
attrition-rate coefficients mentioned above. A comparison is 
made between the approximation in each of these cases and 
the numerical results given by Taylor and Brown. Numerical 
examples and analytic considerations will be utilized to 
investigate the accuracy of the approximation in these 
cases. 



The organization of this thesis is as follows. We first 
review Lanchester’s classic model of combat between two 
homogeneous forces. These results are then extended to the 
variable attrition-rate coefficient case, and the specific 
forms to be studied are discussed. The representation of 
the solution to these variable-coefficient equations in 
terms of so-called general Lanchester functions is 
presented. Then Liouville’s normal form and approximation 
are introduced, evaluated, and discussed for the two 
specific forms of attrition-rate coefficients. 
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II. LANCHESTER'S classic formolation 



F. <J. Lanchester was an English aeronautical engineer 
who first developed a mathematical formulation for combat 
between two conflicting forces. He published this theory in 

a series of articles in the British journal Engineering in 

1914 [Ref 8]. Lanchester's objectives were to provide 
insights into ,, modern ,, combat and to quantitatively justify 
the military principle of concentration. He stressed that 
advances in technology had brought long-range delivery 
capability of modern weapons which allowed for concentration 
of firepower. Lanchester hypothesized that under these 
modern conditions, combat between two military forces could 
be modelled by: 

dx/dt = -ay with x (t=0) = x , 

0 

dy/dt = -bx with y (t=Q) = y^ , (1) 

where x (t) and y (t) are the numbers of X and T forces at 

time t, and t=0 denotes the time at which the battle begins. 

These equations are valid only for x and y greater than 
zero. The parameters a and b are non-negative constants that 
have come to be referred to as Lanchester attrition-rate 
coefficients. These coefficients represent each side's fire 
effectiveness or firepower. The most common set of 
circumstances . under which these equations have been 
hypothesized to apply is that both sides use aimed fire and 
target acquisition times are constant [Ref 16]. 

From (1) Lanchester deduced his classic square law 

2 2 2 2 
b(X Q - x (t)) = a (y - y (t) ) , (2) 

which has the important implication that a force can 
significantly reduce its casualties by initially committing 
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more forces to the battle. The square law also shows that in 
this simple model the Y force wins (ie., X is annihilated) 
if and only if 

x /y < Va/b . ( 3 ) 
0 0 



From (1) Lanchester also derived the now well-known 
results for the time history of the force levels 



x (t) = x coshv/alb t - y x/a/F sinh\/ab t , 
0 0 

y (t) = y cosh Vab t - x x/b/a sinhx/aF t , 
0 0 



(4) 



The victory-prediction condition (3) was given for a 
'•fight to the finish". As H. K. Weiss [Ref 15] points out, 
most battles terminate before complete annihilation of one 
force occurs. Different models of battle termination have 
been proposed, among them that combat ends when either of 
two given "breakpoint" force ratios is reached. Taylor has 
shown that (3) is necessary and sufficient for Y to "win" a 
fixed-force-ratio-breakpoint battle [Ref 14]. 

Lanchester' s original formulation (1) involved many 
implicit assumptions that limit its usefulness for modelling 
real combat. The forces were assumed to be homogeneous, 
there were no replacements or withdrawals, and no movement 
of forces was considered. Setting the attrition-rate 
coefficients a and b in (1) as constants implies that the 
fire effectiveness per firer of each side is constant 
throughout the battle. 



III. VARIABLE AT TRITION-RATE COEFFICIEN TS 

When modelling mobile weapon systems and movement of 
forces, one recognizes that the assumption of constant fire 
effectiveness is open to question [Ref 2]. If fire 
effectiveness is permitted to vary during the battle, 
Lanchester*s formulation (1) takes the form: 

dx/dt = -a(t) y with x (t=0) = x , 

0 

dy/dt = -b(t) x with y (t=0) = y^ , (5) 

where a(t) and b (t) are time-dependent attrition-rate 
coefficients. The form of the attrition-rate coefficients 
depends on such variables as force separation (range of 
battle) , tactical posture of the targets, firing rate and 
rate of target acquisition [Ref 5]. 

New operations research techniques for predicting these 
attrition-rate coefficients were developed in the 1960's. A 
significant contribution in this area was the development of 
a methodology by S. Bonder [Ref 2 and 4] for prediction of 
coefficients based on weapon system performance data. 
Another important development was the method of G. Clark for 
the maximum likelihood estimation of these coefficients 
based on Monte Carlo simulation output data [Ref 6]. The 
work of these and others has facilitated the use of models 
such as (5) in current combat studies. 

To effectively model today's highly mobile warfare, the 
form of the attrition-rate coefficients in (5) must 
accomodate current tactics and battle conditions. 
Specifically, the coefficients must permit the fire 
effectiveness of one or both sides to be non-zero at the 
start of battle (to model, for example, an ambush). The 
coefficients must also allow for the opposing weapon systems 
to have different maximum effective ranges. Since mobile 
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forces close with each other in combat, their fire 
effectiveness should be permitted to increase or decrease as 
a function of range or time. The attrition-rate coefficients 
examined in this paper have the flexibility to model these 
various situations. Motivated by previous work by Bonder 
and Farrell [Sef 5], Taylor and Brown have considered the 
following attrition-rate coefficients: 

m n 

a (t) = k (t+C) and b (t) = k (t+C+A) . (6) 

a b 

The specific cases of (6) to be considered here will be that 
in which A=Q, C>0, 

m n 

a (t) = k (t+C) and b(t) = k (t+C) , m 

a b K ' ' 

and that in which m=n=1, A>0, C>0, 

a (t) = k (t+C) and b (t) = k (t+C+A). (8) 

a b 

The parameter A will be referred to as the '’offset 

parameter” since it permits modelling of combat between 

forces whose weapon systems have different maximum effective 
ranges. The parameter A also allows the ratio of 
attrition-rate coefficients to be non-constant. Koopman [see 
Ref 10] first observed that if the ratio of attrition-rate 
coefficients is constant, the time solution to (5) is no 
more complicated than the constant-coefficient solution. 
That is, when 

a (t) = k h (t) and b (t) = k h (t) , (9) 

a b 

where h (t) is the common time-dependent factor in both 
attrition-rate coefficients, the solution is 

x (t) = x cosh W (t) - y J k /k sinh W (t) , 

0 0 a b 
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where W (t) = \fk k ^lifs) ds. Koopman also found that the 

a b o 

constant ratio of coefficients again yields a square-law 
relationship (see also Ref 16). 

The parameter C in (7) and (8) allows for the initial 
fire effectiveness of each side to be non-zero. This permits 
modelling of battles which begin within the maximum 
effective ranges of the opposing weapon systems. The 
parameters m and n allow the fire effectiveness of each side 
to increase or decrease with time. With m,n<0, the fire 
effectivenesses decrease. This situation might model a 
battle in which combatants take cover and improve their 
positions as the battle progresses. With m,n>0, the 
effectivenesses increase, such as might be the case in which 
two forces are closing with each other. 

To model an engagement in which forces close with one 
another, it is beneficial to relate their instantaneous fire 
effectiveness to the range separating them. As this range 
decreases, one might expect fire effectiveness to increase. 
The attrition-rate coefficients (6) can be expressed in 
terms of range between forces and the maximum effective 
ranges of the opposing weapon systems: 

dx/dt = -a (r) y = -a (1-r/S ) y, 

0 a 

a 

dy/dt - -b (r) x = -b (1-r/R ) x. (10) 

0 b 

These equations have been used by S. Bonder [Refs 1 and 3] 
to model a constant-speed attack on a static defensive 
position. Figures 1 and 2 show sample relationships between 
the attrition-rate coefficients and range between opposing 
forces for various input parameters. In these equations, r 
is the range between the forces, and R and R are the 
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maximum effective ranges of the I and X weapon systems. 



respectively. The constants a 



0 0 



and b represent the fire 



effectiveness of each side at r=0. Range and time are 

related by r(t)=R -vt, where R denotes the opening range of 

0 0 

battle, and v is the constant attack speed. With this 
relationship, the attrition-rate coefficients in (10) can be 
related to (7) and (8) . For (7) we have: 



k = a v/R , k = B v/R , A = (R -R )/v, C = (R -R )/v, 
aOabOB ba aO 



where R > R > R . 
b a 0 



m n 



k = a (v/R ) , k = B (v/R ) , c 

a 0 a b 0 b 




where R = R > R , while for (8) we have: 
a b 0 
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IV. A STANDARD FORM FOR THE VARIABLE-COEFFICIEN T SO LUTION 



The variable attrition-rate coefficient equations (5) 
yield the force level equations 

2 2 

d x/dt - [1/a(t) da/dt ] dx/dt - a(t)b(t) x = 0, (13) 
with initial conditions 



x (t=0) 



x and [ 1/a (t) dx/dt] = -y , 
0 t=0 0 



and similarly 

2 2 

d y/dt - [ 1/b (t) db/dt] dy/dt - a(t)b(t) y = 0 f (14) 

with initial conditions 

y (t=0) = y and [1/b(t) dy/dt] = -x . 

0 t=0 0 

The solution to (13) is given by 

x (t) = C X (t) + C X (t) f (15) 

11 2 2 

where [X (t) , X (t) ] denotes a fundamental system of 
1 2 

solutions to (13) , and C and C are constants determined by 

1 2 

the initial conditions [see, for example. Ref 7]. The 

I-force level has similar form. The functions X (t) , X (t) , 

1 2 

Y (t) , and Y (t) are referred to by Taylor and Brown as 

1 2 

General Lanchester Functions. Their properties are 

summarized in Table I. The functions X (t) and Y (t) are 

1 1 

similar to the hyperbolic cosine, while X (t) and Y (t) are 

2 2 

similar to the hyperbolic sine. The General Lanchester 

Functions may be constructed by successive approximations or 
by infinite series methods. 



16 



A. APPLICATION TO POWER ATTRITION-RATE COEFFICIENTS 



When this solution method is applied to the case of 
power attrition-rate coefficients (7) , the solution (15) may 
be written as 



X(t) s x [0 (C) u (t+C) - V (C)v (t+C) ] - 

0 m,n m,n m,n m,n 

y Jk /k [u (C) v (t+C) -u v (C) ], (16) 

0 a b m,n m,n m,n m,n 

where u , v ,0 and 7 are referred to by Taylor as 
m,n m,n m,n m,n 

Power Lanchester Functions. They have the properties shown 
in Table II. For computational convenience, utilization is 
made of the Lanchester-Clif f ord-Schlaf li (LCS) functions 

oo 2k ^-1 

F (x) = E (x/2) / [k n (j+v) ], 

v k = 0 j = 0 

2k+ 1 k 

G (x) = ? (x/2) / [k n (j + v) ]. 

v k=0 j=0 

The Power Lanchester Functions may be expressed in terms of 
the LCS functions through the following relationships: 



u (t) 
m , n 

v (t) 
m , n 

Q (t) 
m r n 

7 (t) 

m r n 



F (S (t) ) 

q 



(m-n)/2 
t G 

P 



(S (t)) 



F (S(t)) 
P 



(m-n)/2 



G (S (t) ) 

q 



9 



where S (t) 




m n 

(k t ) (k t ) t 
a b 



/ 



( (m+n+2) /2) . 



□sing the LCS 



functions, (16) 



may now be written as 
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x(t) = X fF <J(0))F (J(t)) - ( 1+t/C) (m n>/2 G (J(0))G (J(t)) 

0 p q q p 

- y j/a (t=0)/b (t=0) [ (1+t/C) (m_n)/2 F (J (0) ) G (J(t)) - 

o q p 



G (J(0))F (J(t)) ]. 

p q 



(17) 



B. APPLICATION TO OFPSET LINEAR COEFFICIENTS 



Similar techniques are utilized to apply this general 
theory to the case of the offset linear attrition-rate 
coefficients (8) . Hith these coefficients, a fundamental 
system of solutions to (13) is given by 



(t> =f (t+C) , 



x (t)=g(t+C), y (t) =F (t+C) , y (t) =G (t+C) 
2 1 2 



where 







OO 




d k 






2n 


(2n)! 


n k k 4n-k 


f 


(t) 


= z 


{ 


k 


u /2) 


/ 


} Z B A t 






n= 


0 


a 




b 






k= 0 n 






CO 




d k 






2n+ 


1 


n k k 4 n+2-k 


g 


(t) 


= z 


{ 


k 


V2> 




/ <2n 


+ 1) •' } Z C A t 






n= 


0 


a 




b 






k=0 n 






OO 




d k 






2n 




n k k 4n-k 


F 


(t) 


= z 


£ 


k 


u /2) 


/ 


(2n) ! 


} ZD A t 






n= 


0 


a 




b 






k=0 n 






OO 










2n+ 


1 


n+ Jk k 4n+2-k 


G 


(t) 


= z 


c 


(^k 


k 


, /2) 




/ (2n 


+ 1) ! } Z E A t 






n= 


0 


a 




b 






k=0 n 



The coefficients 



k 

B , 
n 



k 

C , 
n 



k 

D , 
n 



k 

and E satisfy complicated 
n 



recurrence relations [see Taylor and Brown (Ref 13) ]. The 

functions f (t) , g (t) , F (t) , and G (t) are referred to as 
Offset Linear Lanchester Functions. As in the case of the 
power attrition-rate coefficients, use is made of the 
following Auxiliary Offset Linear Lanchester Functions for 
computational convenience: 
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h (p r u) 



oo 2n n k k 

= MP / (2n) !} E B p 
n=0 k=0 n 

oo 2n+ 1 n k 1c 

tf(p r u) = E CP / (2n+1) i } SCu 

n=0 k=O n 

oo 2n n k k 

H(p,u) = Z {p / (2n) ! } Z D u 
n=0 k=0 n 

co 2n+ 1 n+1 k k 

W (P,u) = Z (P /(2n+1) ! } Z E u . 

n=0 k=0 n 

Otilizing these auxiliary Lanchester Functions, the X-force 
level solution (15) may be written as 

x(t) = x [H (J(0) ,j (0) ) h (J (t) , j (t) ) - W (J (0) (0) ) 

0 

j (t) ) ] 

- yjk /k ( h (J (0) ,j (0)) w (J (t) , j (t) ) - w(J(0) ,j(0)) 

0 a b 

h(J(t) ,j(t)) ], (18) 

f 2 

where j (t) = A/(t+C) and J(t) =\l k k (t+C) /2. 

a b 

The equations (17) and (18) are solutions to the 
X-force level equation (13) for the case of power and offset 
linear attrition- rate coefficients, respectively. This 
solution method of Taylor and Brown represents a 

considerable extension of Lanchester theory in the area of 
variable a ttrition-rate coefficients. Previous analytic 
results had been available only in the more limited case of 
the opening range of battle equal to the minimum cf the 
maximum effective ranges of the two weapon systems. That is, 
the initial fire effectiveness of at least one side was 
restricted to be zero [see for example Bef 12]. 

While the time solutions (17) and (IB) appear 

formidable, their implementation is simplified if a digital 
computer is available. As noted earlier in this paper. 
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however, these solutions are complex and involve new 
mathematical functions. This makes it difficult for the 
military analyst to perceive significant relationships among 
battle variables, and renders parametric analysis difficult 
and time-consuming. This provides motivation for the 
approximation technique to follow. 
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V. 



REDUCTION TO LIOUVILLE^ NORMAL FORM AND APPROXIMATION 



The argument ab (referred to by, among others, Taylor 
and Brown [Ref 13] as the "intensity of combat") in the 
constant coefficient result (4) provides motivation for the 
transformation [Ref 14] 

t 

z = f\la. (s) b (s) ds, (19) 

c 0 

where z (t=0) is denoted z . Applying this transformation to 

0 

(13) and (14) yields 
2 2 

d x/dz + {1/2 d/dz ln[b(t)/a(t) ]} dx/dz - x = 0, (20) 
with initial conditions 



x ( z=z o ) 



1/2 

x and {[ b (t) /a (t) ] dx/dz} = -y , 
0 z=z 0 



and 



2 2 

d y/dz * {1/2 d/dz ln[a(t)/b(t) ]} dy/dz - y = 0, (21) 
with initial conditions 

1/2 

y (z=z ) = y and {[a(t)/b(t)] dy/dz} = -x . 

0 0 z=z 0 



Taylor and Brown have shown that the X force level 
equation (13) may be transformed into a linear second order 
differential equation with constant coefficients if and only 
if (20) is a constant coefficient equation. This is the case 
if and only if 

(1/Va (t) b (t) ) d/dt ln[a(t)/b (t) ] = CONSTANT. 

Equation (20) shows the significance of the parameters 
a(t)/b(t) (the relative effectiveness of the two weapon 
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systems) , and the intensity of combat \la. (t) b (t) in 
determining battle outcome. This is a generalization of the 
well-known constant coefficient results to the case of 
temporal variations in fire effectiveness. The significance 
of these parameters may be seen more explictly by 
transformation to Liouville's normal form (see Ref 14). 



The transformation 

1/4 

x(z) = X(z)[ (a(t)/a ) / (b(t)/b Q )] , (22) 

may be applied to (20) and yields the normal form with the 
first derivative of the dependent variable removed 



2 2 

d X/dz - [ 1 + F (z) ]X = 0, 



(23) 



with initial conditions 



X (z=z ) = x and dx/dz(z=z ) = -y ya. /b - x e , 
00 0 00000 



where 



i • 



P(z) = [ R (t) ] 



-1/4 



(24) 



F(Z) = P (Z) / P (Z) 

R (t) = a (t) /b (t) e (t) = 1/4 \/a(t) b (t) d/dt InR (t) . 

(25) 

2 2 

In (24), the notation P''(z) denotes d P/dz . Equation (23) 
is Liouville's normal form (see page 23 of Ref 9]. 



For the Y force level, the transformation 

1/4 

y(z) * I(z)[ (b(t)/b ) / (a(t)/a Q )] 
yields the normal form 



2 2 

d Y/dz - [1 + G (z) ]Y = 0, 



(26) 



with initial conditions 
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Y (z=z ) = y and dY/dz(z=z ) = -x yb/a + y e , 

00 0 00000 



where 



1/4 

G (z) = Q" (z) / Q(z) and Q(z) = [ R (t) ] = 1/P (z) 



Expressing the transformed force-level equation (23) in 
the simple form 



2 2 

d X/dz - X = Q (z) X, 

variation of parameters may be used to obtain the solution 
to (23) as 

X (z) = x cosh (z-z ) - [y \/a /b + x e ] sinh (z-z ) 

0 0 00000 0 

+ f? ( j) sinh (z- j) X ( j) d j . (27) 

z o 

This equation may be expressed in terms of the original time 
variable t and dependent variable x as 

1/4 t l 

x(t) = [ (a(t)/a )/(b(t)/b ] {x cosh ( /& (s) b (s) ds) 

0 0 t 0 0 

-[ y \J a /b + x e ]sinh ( a (s) b (s) ds) } 

0 0 0 0 0 o 

1/4 1/4 

+ C a (t) /b (t) ] d/£s{1/(a(s)b(s) ) d/ds[ b (s) /a (s) ] } 

sinh ( /i/a (sj) b ( j ) dj) X (s) ds. (28) 

s 

If the relative effectiveness varies slowly, then the 
relationship (24) would indicate that F(z)<<1 so that one 
could drop the integral term in (27) yielding 

X (z) = x cosh (z-z ) - [y \Ta/b + x e ] sinh (z-z ), 

0 0 00000 0 { 29 ) 

which is the Liouville-Green approximation [Ref 11]. The 
approximation may be expressed in terms of the original 
variables t and x as 
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X(t) 



1/4 t 

[ (a (t)/a )/b (t)/b ] [x cosh(/Va(s) b (s) ds) 

0 o 1 0 0 

- [7 >/a /b + x e ]sinh ( /^/a (s) b (s) ds) } (30) 

0 0 0 0 0 o 
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VI 



APPLICATION TO POWER AND OFFSET LINEAR CO EFFICIENTS 



The Liouville-Green ap proximation may now be applied to 
the constant-speed attack model with attrition- rate 
coefficients (7). For the reader's convenience, these 
coefficients are 



m 

a (t) = k (t+C) 
a 



and 



b (t) =? k (t+C) Q . 

b 



(7) 



This model assumes that the maximum effective ranges of the 
two opposing weapon systems are equal (no offset) . The 
opening range of battle is taken to be within this maximum 
effective range (i.e., C>0) . Applying the Liouville-Green 
approximation (29) , the time solution to the X force level 
equation in terms of the original time variable t and 
dependent variable X becomes 



(m-n)/4 c . 

x(t) = (1+t/c) (x cosh/va (s) b (s) ds 

0 0 

. (m-n) /2 -(m+n+2)/2 

- [y Vk /k C + (x (m-n)/4 k k )C ] 

0 a b t 0 a b 

sinh jVa (s) b (s) ds) , 

0 



(31) 



where 



/\/a (s) b (s) ds = (V k k /(m+n + 2)/2) 
0 a b 

(m+n+2) /2 (m+n+2)/2 

{ (t+C) -C } . 



For future reference, equation (24) becomes in this 

case 

2 2 

F (z) = [ (m-n) (3 m + n+4) ] / [4(m + n + 2) z ], (32) 

or in terms of the original time variable t and range 
parameters (11): 
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m n 

F(t) = [fi R (m-n) (3m + n+4) ] / 
a b 

m + n m+n+2 

[16a b (vel) (t+(R -R )/vel) ], (33) 

0 0 a 0 



where vel = constant velocity of attack. 



We can also apply the Liouville-Green approximation to 
the case of the offset linear attrition-rate coefficients 



a (t) = k (t+C) and b (t) = k (t+C+A) . (8) 

a b 



In this case, the time solution to the X force level 
equation beccmes 



1/4 } 

x(t) = [ (1+ (A/C)/(1+ (A/ (t+C) ) ) ] {x cosh f\]3L (s) b (s) ds 

o 0 

- [ (y \/k /k )/(l+ (A/C)) + (X (A/C))/ 

0 a b 0 

2 3/2 c 

(4 k k C ( 1+ (A/C) ) ) ]sinh/Va (s) b (s) ds} , (34) 

a b 0 



where 



/Va (s) b (s) ds = (A k k /8) { (1+ (2 (t+C) /A) ) (1+ (2 (t +C) /A) ) -1 

0 a b 

2 

- ln[ 1+ (2 (t+C) /A) + (1+ (2 (t+C) /A) ) -1]} 



For future reference, equation (24) becomes in this 

case 



F (t) = R R (R -R ) [ 12 (t + (R -R ) /vel) + 7 ( (R -R )/vel)] / 
aBBa aO Ba 

3 3 3 , % 

{16a b vel (t+ (R -R )/vel) (t+ (R -R )/vel) }. (35) 

0 0 a 0 b 0 
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VII. DISCUSSION 



A computer program was developed in FORTRAN to produce 
plots of the X force level versus range between opposing 
forces for the constant-speed attack model described 
earlier. Both the power and offset linear attrition-rate 
coefficient cases (equations (31) and (34)) were considered. 
These approximate solutions were compared to the exact 
results given by Taylor and Brown [Ref 13] to evaluate the 
accuracy of the Liouville-Green approximation to equations 
with these coefficients. 



A. POWER ATTRITION-RATE COEFFICIENTS 

For the case of power attrition-rate coefficients, the 

maximum effective ranges of the two weapon systems were set 

equal to 2000 meters, with an opening range of battle of 

1250 meters. The fire effectiveness of the two sides at zero 

range separation was a =.06 X casualties/min-Y unit, and 

0 

b =.6 Y casual ties/min-X unit. The velocity of attack was 
0 

constant at five miles per hour, and the initial force 
levels were X = 10 and Y = 30. These parameter values 

correspond to those used by Taylor and Brown in their paper 
[Ref 13]. The exponents m and n in (7) were varied 
throughout the range explored in Ref 13. The results are 
shown in Figures 3 through 10, where each figure contains 
the approximate time solution and the corresponding exact 
analytical solution obtained by Taylor and Brown. As even a 
cursory examination reveals, the approximation appears to 
conform to the established results only in certain cases. 
Some possible explanations for this observed deviation will 
now be considered. 

Recall that the integral term in equation (27) was 
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dropped (assumed to be approximately zero) in making the 
Liouville-Green approximation. The validity of this 
assumption appears to be linked strongly to the form of 
F (z) , which for the present case is 

m n 

F ( t) = [R R (m-n) (3m + n+4) ] / 
a b 

m+n m+n +2 

[16a Q b o (vel) <t+ (R^-R^ /vel) ]. ( 33 ) 

When the exponents m and n are equal, F(t)=0, and by 
considering equation (27) , one can see that the 
approximation is, in fact, the exact solution. This is 
verified trivially by Figures 3 and 4, for which the 
exponents are equal. It should be noted that in this 
situation, the model reduces to the case of a constant ratio 
of attrition-rate coefficients. 



When the exponent m is greater than n, F (t) will always 
be greater than zero, and from (27) it follows that the 
exact solution is greater than the approximation. This 
situation is illustrated by Figures 6 and 10. The converse 
also holds. When m is less than n, F(t) is always less than 
zero, and the approximation is greater than the exact 
solution. This is shown in Figures 5 and 9. 



Again examining (33) 



the parameter a b , which is the 
1 0 0 1 



product of the opposing forces* fire effectivenesses at r=0, 
appears to be significant. As |a b | increases, F(t) 



decreases, making the integral term in (27) smaller, so one 



would expect the approximation to be "closer" to the exact 
solution. This did occur in the cases considered (see for 
example Figure 7) . 



B. OFFSET LINEAR COEFFICIENTS 



28 



For the case of the offset linear attrition-rate 

coefficients (8) , the maximum effective range of the Y force 

weapon system (R ) was set equal to 1500 meters, with the 

a 

rest of the parameter values equal to those in the power 

coefficient case. The maximum effective range of the X force 

weapon system (R ) was varied throughout the range 

b 

considered in Ref 13. The results are shown in Figures 11 

through 15 where, as before, each figure compares the 
corresponding solutions.. As was the case with the power 
coefficients, the solutions are seen to conform only in 
certain cases. 

For the linear coefficients, it was shown that F (z) 
could be written in terms of the time variable t as 

F(t) = R R (R -R )[12(t+(R -R ) /vel) + 7 ( (R -R )/vel) ] / 
a3Ba aO Ba 

3 3 3 , x 

{16a b vel (t+ (R -R ) /vel) (t+ (R -R )/vel) ). (35) 

0 0 a 0 b 0 

When the maximum effective ranges of the opposing weapon 

systems (R and R ) are equal, F (t) =0, and the integral 
a b 

term in (27) is zero. This again reduces to the constant 

attrition-rate-ratio case as can be seen from (8) with A=0. 
Thus one would again expect agreement of the exact and 
approximate solutions. This was found to be true as shown 
by Figure 11, where both maximum effective ranges equal 1500 
meters. 

As was the case for the power coefficients, the 
parameter |a^b^| appears in the denominator of F (t) , hence 

the two solutions should be closer as la b increases. This 

1 0 0 1 
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was found to be the case, as illustrated by Figure 15. 



C. ANALYTICAL CONSIDERATIONS 

The examinations of both the power and offset linear 

attrition’- rate coefficient cases have centered on the 

parameter F (z) in the solution (27). While F (z) does appear 

to be the driving factor in the integral term in (27) (thus 

in the validity of the approximation) , it should be made 

clear that the other factors in the integral term can not be 

ignored. While a mathematical analysis (for example of the 

hyperbolic sine term) would be formidible, examples can be 

found which illustrate that F (z) can not always be 

considered alone. One such example is shown in Figure 8. The 

fire effectiveness of the two forces at r=0 were 

interchanged (i.e., a =.6 and b =.06 instead of vice versa) 

0 0 

in the power attrition-rate coefficient case. From the 

appropriate form of F (z) , equation (33), one can see that 
this change has no effect on the numerical value of F(t) 
throughout the time range. Yet the significant difference in 
results can be seen in Figure 8 when compared to Figure 6. 
Currently, there is no convenient methodology to analyze the 
implications of such discrepancies. 

Toward this end, theoretical analysis in the area of 
Olver's “Error Bounds for the Liouville-Green Approximation” 
[Ref 11] could be very beneficial. While they sound 
promising, Olver's results can not be directly applied to 
variable-coefficient Lanchester-type differential equations. 
If available in a convenient form, such error bounds could 
be evaluated for specific situations to investigate 
tolerable error limits prior to invoking the approximation. 

From the above analytic considerations, one can see 
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that while one may be able to explain some of the error 
incurred, this is small comfort to the analyst attempting to 
make practical use of the Liouville-Green approximation. It 
has been shown that the approximation appears to yield 
credible results only in certain limited circumstances. 
Furthermore, methodology is not presently available to 
predict the error which will be incurred in a specific 
situation, except for the trivial case in which the integral 
term in (27) goes to zero. 
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VIII. 



SOHMABY 



This thesis has examined the "adequacy" of the 
so-called Liouville-Green approximation to the solution of 
variable-coefficient Lanchester-type equations of modern 
warfare. Although this approximation only involves 
"elementary" functions and has an intuitively appealing 
form, it was unfortunately found that this 
Liouville-Green-Lanchester approximation is not consistently 
reliable for # estimating force levels. Specif ically, the 
approximation was applied to two forms of attrition- rate 
coefficients: (I) power attrition-rate coefficients 
(modelling, for example, the same maximum effective ranges 
for the opposing weapon systems) , and (II) offset linear 
attrition-rate coefficients (modelling different maximum 
effective ranges) . It was seen that the approximation was 
made by dropping the integral term in equation (27) , hence 
exact results are achieved only when this term equals zero. 
This situation was found to occur only when the ratio of 
attrition-rate coefficients was constant. For the power 
attrition-rate coeff icients, this implies that the exponents 
in (7) are equal, while for the offset linear coefficients 
this implies the offset parameter A in (8) is zero (which 
occurs when the maximum effective ranges are equal) . 
Analytic considerations showed that while certain initial 
parameter combinations produced credible results from the 
approximation, there was no methodology available to predict 
the error incurred in the general case. It was found that, 
as a general rule, the further one moves from the constant 
coefficient ratio case, the larger the error incurred by 
invoking the Liouville-Green approximation. 
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TABLE I. 



Properties of the General Lanchester Functions x^, x^j y^, y 2 - 

1. dx^/dt = /k^/k^ a(t)y 2 



2. dx 2 /dt = A ~/T a(t)y x 



3. dy 1 /dt = /k a /k b b(t)x 2 

4. dy 2 /dt = * /l k^7k^ b(t)x 1 

5. J x 1 (t)y 1 (t) - x 2 (t)y 2 (t) =1 Vt 
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Table II. Properties of the Power Lanchester Functions 



u , v , U , and V 
m,n m,n m,n n 



1 . du /dt = ^k. /k [k t m ]V (t) 

m,n b a a rn,n 



2 . 



dv /dt = /k, /k [k t m ]U (t) 
m,n d a a m,n 



3 * dU /dt = A/k. [k, t n ]v ft) 

m,n a b b 01,11 

4 . dV /dt = /k ~ 7 kT • [ 1 c t n ]u (t) 

m , n a b b 01,11 



5. 



u _(t)u (t) - v (t)V (t) 

tn,n 01, n m, n 01, n 



6 . 



u (t= 0 ) = U (t= 0 ) = 1 
01, n tn,n 



7. 



v (t= 0 ) = V (t= 0 ) = 0 
m,n m,n 



8 . u (t) = U (t) = cosh (/k k, 

m,m 01,111 a b 



9. V (t) = V (t) = sinh (/k k, 

01,01 m,m a b 



Properties 6 and 7 only hold for m,n > - 1 . 



= 1 Vt 

t m+ l/ (nH-1) ) 
^/(nri-l)) 
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of battle Is denoted as R~ •» 1250 meters and (as shown) 



